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We investigate the connection between the NSVZ and the DRED forms of the gauge 
/^-function in an iV = 1 supersymmetric gauge theory. We construct a coupling constant 
redefinition that relates the two forms up to four loops. By abelian calculations, we are 
able to infer the complete non-abelian form of [3^ DRED , and also (3g 4 ^ DRED except for 
one undetermined parameter. 
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1. Introduction 



An all-orders formula for the gauge /5-function in an iV = 1 supersymmetric gauge 
theory was presented some years ago. This result (which we shall call (3^ svz ) originally 
appeared (for the special case of no chiral superfields) in Ref. |J, and was subsequently 
generalised, using instanton calculus, in Ref. 0. (See also Ref. 0.) For a recent discussion 
emphasising the importance of holomorphy, see . 

Recently the renormalisation group fixed points of j3^ svz have become important 
in the study of duality (for a review, see Ref. An interesting question, therefore, is 
as follows: given the renormalisation scheme dependence of /3-functions beyond one loop, 
in which scheme is the NSVZ result valid? For instance, will calculations using standard 
dimensional reduction (DRED) give the NSVZ result? The DRED result certainly agrees 
with the NSVZ result at one and two loops; moreover certain properties of (3 g DRED at 
higher loops are consistent with the NSVZ result. Namely, [3® RED is known to vanish 
at three loops for a one-loop finite theory || J7| and furthermore, if we specialise to the 
N = 2 case, (3 g and the anomalous dimensions of the chiral superfields vanish beyond one 
loop ||. One might accordingly be tempted to speculate that DRED will reproduce the 
NSVZ formula to all orders. However, in a recent note |J we showed that this is not the 
case; at three loops the DRED result is related to the NSVZ result by a coupling constant 
redefinition. In the present paper we shall give more details of this calculation and also 
extend the result to four loops, at least in the abelian case. 

Before proceeding, however, it is worthwhile emphasising the following point. It is 
sometimes asserted that the perturbative coefficients of (3 g are quite arbitrary beyond two 
loops, so that, for example, all contributions at three and more loops can be transformed to 
zero. We shall see, however, that in the general case (with a superpotential) the nature of 
possible changes in (3 g due to redefinitions 5g which are manifestly gauge-invariant analytic 
functions of g and the Yukawa couplings Y l ^ k is heavily constrained. (One has also the 
freedom to make redefinitions 5Y l i k , but, as we shall see, these are not germane to the 
issue of whether the NSVZ and DRED results are equivalent.) 

1 We will discuss later whether this result will apply in schemes other than DRED. 
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2. The three-loop calculation 

The Lagrangian -£/susy(W) for an TV = 1 supersymmetric theory is defined by the 
superpotential 

W = \Y iik ^^ k + (2.1) 

-^susy is the Lagrangian for the N = 1 supersymmetric gauge theory, containing the 
gauge multiplet V and a multiplet of chiral superfields $^ with component fields {4>i,ipi}, 
transforming as a representation R of the gauge group Q. We assume that there are no 
gauge-singlet fields. The /3-functions for the Yukawa couplings 0^ k are given by 

0Vk = Y P^ 1 k ) p = Y ijp 7 k p + (k^i) + (k^ j), (2.2) 

where 7 is the anomalous dimension for $. The one- loop results for the gauge coupling 
/3-function j3 g and for 7 are given by 

\QTx 2 [5 g l) = a 3 Q, and 167r 2 7 (1)i j = P^, (2.3) 

where 

Q = T(R) - 3C(G), and (2.4a) 

pi. = i Y ikl Y jkl - 2g 2 C(RY J . (2.46) 

Here Y jkl = (Y^ kl )*, and 

T(R)5 AB =tv{R A R B ), C(G)5 AB = f ACD fBCD and C(RY J = (R A R A y j . (2.5) 

The two-loop /3-functions for the dimensionless couplings were calculated in 
Refs. 101, PHl3fl: 



(16tt 2 ) 2 ^ 2) = 2g 5 C(G)Q - 2g 3 r - 1 C(RY 3 P^ (2.6a) 
WjV 2 ^ = [-Y jmn Y mpi - 2g 2 C{RY 3 5\]P n p + 2g*C(RY J Q, (2.66) 



where Q and P l j are given by Eq. (|2.4|) , and r = S AA . 
In our notation the NSVZ formula for j3 g is 



a 3 

nNSVZ __ il 
Pg 16tT S 



Q-2r-Hr [jC(R)} 
l-2C(G)g 2 (16n 2 )~ 1 



(2.7) 
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which leads to 

{l^ 2 f(3^ NSVZ = 4g 7 QC(G) 2 - 4^ 5 C(G r )r -_1 167r 2 tr 1^0(11) 
-2f7 3 r- 1 (16vr 2 ) 2 tr W 2) C(R) 



(2- 



As mentioned in the Introduction, one would like to know in which renormalisation scheme 
this result is valid, since /3-functions are scheme dependent. It is easily seen from Eqs. ( |2.3| ) 
and ( |2.6| ) that the NSVZ result coincides with the results of DRED up to two loops. Will 
DRED reproduce the NSVZ result to all orders? There are two pieces of evidence which 
appear to favour this conjecture. Firstly, ^^ DRED h as been explicitly shown to vanish 
for a one-loop finite theory ||, i.e. one for which P = Q = 0, and it is clear from Eqs. ( |2.8| ) 
and (|2.6| ) that this property is shared by [5^ NSVZ . A second piece of evidence comes from 



specialising to the N = 2 case. In N = 1 language, an N = 2 theory is defined by the 
superpotential 

W = VtgriAtfSAUZj ( 2 -9) 

where 77, x an d £ transform according to the adjoint, S* and S representations respectively. 
The set of chiral superfields x, £ is called a hypermultiplet. N = 2 theories have one-loop 
divergences only||]; using DRED we may therefore expect that (3 g and the anomalous 
dimension of both the r\ and the hypermultiplet should vanish beyond one loop. I We see 
that the NSVZ result of Eq. ( |2.7|) is consistent with this property; we have 

PriAB = Qg 2 ^AB 

(2.10) 

so that if 7 vanishes beyond one loop then Eq. ( |2.7| ) reduces to (3^ svz = j§^zQ, which 
is of course the one-loop result. In particular, since from Eq. (|2.66| ) we have 7^) = for 
iV = 2, we have (5^ NSVZ = for an iV = 2 theory. Nevertheless, despite these indications 
that the NSVZ result might coincide exactly with that obtained using DRED, we shall 
now show that in fact the NSVZ and DRED results part company at three loops. We shall 
see that they are related by a coupling constant redefinition corresponding to a change of 
renormalisation scheme. We shall calculate ^^> DRED explicitly in the abelian case and 



2 Clearly the absence of divergences beyond one loop implies that f3 g vanishes beyond one loop 
as long as minimal subtraction is employed; higher order contributions to the (3 g can be invoked 
by making finite subtractions, or equivalently by a redefinition g — >• g + 8g. 
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construct the coupling constant redefinition which effects the transition to the NSVZ result. 
We shall then extend this redefinition to the full non-abelian case by exploiting the known 
N = 2 properties of the DRED result, and hence we shall deduce the full non-abelian 

n(3)DRED 
Pg 

We calculate (3 g by computing the divergences in the vector field two-point function, 
using the super- Feynman gauge. This is sufficient since in the abelian case the background 
superfield calculation and the normal superfield calculation are identical. The relevant 
diagrams are shown in Fig. 1. The shaded blobs represent one-loop self-energy insertions. 
We use dimensional reduction (DRED) with minimal subtraction, setting e = 4 — d. The 
divergent part of each individual diagram, after performing all the appropriate subtractions 
for divergent sub-diagrams, will be expressible as a combination ATli + _Bllo, where A and 
B are analytic in - and the coupling constants, and 111 and Ilo are projection operators 
defined in the Appendix, with 111 + Ho = —d 2 . Upon adding the results for all diagrams, 
the total will only involve IIi, reflecting the transversality of the vector propagator. We 
use a convenient and efficient short-cut to calculate the total coefficient for II i without 

2 

calculating the full contribution for each diagram. The idea is the following: To obtain the 
total coefficient of II i , it will be sufficient to know the difference B — A for each diagram. 
Upon summing over all diagrams, the sum over the Bs will give zero just leaving the sum 
over the As, which is what we want. The point is that we can obtain the combination 
2(5 — A) simply by adding up the divergent contribution to D 2 D 2 from each diagram, 
regarding the derivatives as exactly anticommuting. (So that D 2 D 2 , D 2 D 2 and D a D 2 D a 
would each count the same.) Each diagram may start with up to 10 Ds and 10 Ds (all 
on internal lines). For each diagram, we manipulate the supercovariant derivatives D 
and D using integration by parts and the anticommutation relations Eq. (A.l) until we 
obtain a set of diagrams each of which contains 8 -Ds and 8 -Ds, possibly together with 
some ordinary derivatives. During this process we avoid integrating any D or D onto an 
external line. The 3 cfiO integrals will absorb 6 .Ds and 6 .Ds, leaving 2 Ds and 2 Ds 
which will contribute to II i or Ilo. We are only interested in knowing the contribution 
to D 2 D 2 , so from this point on we can treat the Ds and Ds as exactly anticommuting. 
We now use integration by parts to arrange for each loop to contain 2 Ds and 2 Ds, and 
integrate the remaining 2 .Ds and 2 Ds onto an external line, writing them in the form 
D 2 D 2 . We can now do the 9 integrals, leaving us with a momentum integral. We evaluate 
the momentum integral, subtract its subdivergences using minimal subtraction and finally 
obtain the divergent contribution to D 2 D 2 . The reader might like to check the extent to 
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which this trick simplifies the one-loop (or even two-loop) calculation. The disadvantage, 
obviously, is that we lose the nice check afforded by the cancellation of the Ilo terms. 

The results calculated according to the above procedure for each diagram in Fig. 1 
are shown in Table 1. The momentum integrals for each diagram can be expressed (using 
integration by parts) in terms of a basic set depicted in Fig. 2, which can be evaluated 
with their subtractions once and for all. In Fig. 2, a dot on a line represents a squared 
propagator, and arrows represent linear momentum factors in the numerator, the momenta 
which correspond to a pair of arrows being contracted together. Each diagram in Fig. 2 
represents a Feynman integral of logarithmic overall divergence; there are no external lines 
because we have set the external momentum zero. This means that at least one propagator 
must be given a mass, since otherwise there is no scale defined and of course dimensional 
regularisation cannot be employed: the overall logarithmic ultra-violet and infra-red di- 
vergences would cancel. (Alternatively an external momentum may be "threaded" in an 
arbitrary way through the diagram). Also, any explicit infra-red divergence corresponding 
to a squared propagator must also be regularised by introducing a mass; such a mass will 
often serve to define the scale, too. Sometimes an alternative 5-function infra-red regulator 
can be useful [|14| , whereby instead of 

' ' or * (2.11) 



(k 2 ) 2 (k 2 + m 2 ) 2 k 2 (k 2 + m 2 ) 

one has 



' ' +T 1 ~J W (k). (2.12) 



(k 2 ) 2 (k 2 ) 2 4-d 

(Note that we have chosen to perform the diagrammatic calculation in Euclidean space, 
for which bookkeeping of factors of i etc. is easier.) After subtraction of ultra-violet 
sub divergences, the result for each diagram is independent both of the means by which a 
scale is introduced and the method used to regulate infra-red divergences. (Note that this 
statement is true only of the subtracted diagram.) The second column of Table 1 shows the 
expression for each momentum integral in terms of those in Fig. 2. (The presence of a zero 
in this column may indicate that the -D-algebra for the diagram gave zero; or alternatively, 
that the momentum integral reduced to a "factorised" form which can be shown JT5| to 
produce no simple pole after subtraction.) Each diagram in Fig. 1 also corresponds to 
a product of group matrices and Yukawa couplings which may be expressed in terms of 
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basic invariants, and which are given in the third column of Table 1 (including also the 
symmetry factor for the diagram). The invariants Xi, X2, X3 and X4 are given by 

Xx = g ^Y klm P n l C{RY m Y knp = ghi [S 4 C(R)} 

X 2 = g 4 Y klm C(R) n iC(R) p m Y knp = gHr [S^R)] , (2.13) 
X 3 = gHr[PC(R) 2 }, X 4 = ghr[P 2 C(R)}. 

(See Appendix A for the definitions of Si, S4 and other notational details). Upon adding 
all the results, we find 

< VV > po ie= 2r~ l {{2A - C){X 1 + 2X 3 - 2g 6 Qtr[C(R) 2 }) + (2B + C - 2D)X A ) . 

(2.14) 

We note that all contributions involving the diagram E have cancelled; this will be impor- 
tant, since E is the only diagram out of our basic set in Fig. 2 whose simple pole involves 
C(3). We also note that the invariant X2 does not appear in the final result, and therefore 
< VV > P oie (and hence j3 g ) vanish in the one-loop finite case, when P = Q = 0.0 It 
is interesting that we can draw this conclusion without evaluating any of the Feynman 
integrals. The /3-function is derived from the simple poles in the subtracted diagrams. We 
find 



A - 4 1 b - ^ - 

^simple 3(167r2) 3 e ' Simple 3 (^2)3,' 

(2.15) 



O simp le - g (167r2)3e > Simple - 3 ^2^3^ 

and hence we have (note that with our conventions, at L loops f3g L ^ differs from the 
corresponding simple pole contribution to < VV > po i e by a factor of ^-) 

(16n 2 ) 3 (3^ DRED = r~ x g {3Xi + 6X3 + X A - 6g G Qtv[C(R) 2 ]} , (2.16) 

(3) 

where r = 5aa- On the other hand, the NSVZ result for f3 g in the abelian case, obtained 
by setting C(G) = in Eq. @, is simply pf )NSVZ = -2 T f^r- 1 tr[ 7 ( 2 )C( J R)] , and so 
we have 

(leTr 2 ) 3 ^ 3 ^ 5 ^^ = r- x g {2X X + 4X 3 - 4g 6 Qtr[C (R) 2 } } . (2.17) 

We note that X4 appears only in the DRED result Eq. ( |2.16| ) , while the other terms appear 
in the same combination X x + 2X 3 - 2g 6 Qtr[C(R) 2 } in both Eqs. (pT6Q and (|2T7|) . We 



3 Of course there is no non-trivial one-loop finite abelian theory; we are here really speaking 
of the corresponding terms in the full non-abelian case. 
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also note that the same combination already appeared in Eq. ( 2.14 ) before evaluating the 
momentum integrals. 

We now wish to show that the two results Eqs. (|2.16| ) and ( |2.17|) are in fact related 
by a coupling constant redefinition, equivalent to a change of renormalisation scheme. In 
general, a redefinition Sg of g induces a change S/3 g of f3 g given (to lowest order in g) by 

SP 9 {g,Y) 



d d d 

dY +(3 * Y -dY^ +P9 -dg~ 



5g-5g.j-[3 g . (2.18) 



In particular, if we choose 



Sg = -(167r 2 )- 2 ir-ytr [PC{R)\ , (2.19) 
then the resulting 8(3 g is given by 

(lGvr 2 ) 3 ^ = r~ l g (-X x - 2X 3 - X 4 + 2g 6 Qtr[C(R) 2 ]) . (2.20) 
We easily see that then 

P&NSVZ = 0&DRBD + gp gm (221) 

It is obvious from this analysis that, as we already mentioned in the Introduction, it is 
quite non-trivial that /3 g 3 ^ NSVZ and p g 3 ^ DRED ma y b e related in this way. Hence this 
provides a strong check on the validity of the NSVZ result. 

In these considerations we explored redefinitions of the gauge coupling g only; what 
of redefinitions of the Yukawa couplings Y 1 ^? Under an arbitrary redefinition 

yijk _^ fyijky (2.22) 

we have at once that 

P g '(Y',g)=P g (Y,g) and j'(Y',g) = 7 (Y i9 ). (2.23) 

It follows that if P g (Y,g) and j(Y,g) satisfy Eq. fl277|), then P g '(Y,g) and j'(Y,g) do 
likewise. This means that a redefinition Y —>■ Y' has no effect on the question of whether 
(P g , 7) obey the NSVZ condition. This does not mean that redefinitions of Y have no 
significance, however; for example, in Ref. [IB], we showed how to construct a redefinition 



of Y corresponding to a change to a scheme such that 7^ 3 - ) = for a one-loop finite theory. 
We will return to this redefinition in the next section. 
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We turn now to the non-abelian case. The crucial observation is that Sg as defined in 
Eq. ( p.!9| ) does not vanish for an N = 2 theory in general (though it does in the abelian 
case, as may be easily verified). There is, however, an obvious generalisation of it to the 
non-abelian case, to wit 

Sg = (IGtt 2 )- 2 ^ 3 [r-hi [PC(R)\ - g 2 QC(G)] (2.24) 

where we have reversed the overall sign (compared to Eq. (|2.19|) ) because we plan to use 
this Sg to go back from f3 g ^ NSVZ to f3 g ^ DRED . It is easy to verify that Eq. ( |2.24j ) leads 
to Sg — in the N = 2 case. Is this the only possible extension of Sg to the non-abelian 
case? We are constrained by the following requirements: 

(1) Sg = for a one-loop finite theory. This is because we know that both f3 g 3 ^ NSVZ 
and p^ DRED vanish in the one-loop finite case. This is manifest in the NSVZ case, 
and was explicitly verified in Ref. [|J for the DRED case. (See also Ref. 0.) In the 
three-loop case the relevant /5-functions in Eq. (|2.18[ ) are one-loop. Since these vanish 



in this case, to produce a vanishing S/3 g , Sg must vanish also. 
(2) Sg = for a iV = 2 theory. In the N = 2 case we know that, as discussed earlier, 
p(3)NSVZ _ p(3)DRED _ g Q ear i v we therefore require S(3 g = 0; Sg = ensures this, 



Py-qy +Py-8y* ^fig-do 



just corresponds to multiplication by 



since the operator 
a factor in the N = 2 case. 

(3) Eq. (|2.19|) must hold in the abelian case. 

(4) The resulting terms in 8/3 g must correspond to possible 1PI Feynman graphs. 

It is easy to convince oneself that Eq. ( |2.24j ) represents the only possible transformation 
(up to an overall constant, which we have fixed by the abelian calculation). Our result for 
non-abelian case is therefore H: 



(167v 2 ) s p g 3)DRED = r~ 1 g{3X 1 + 6X 3 + X 4 - 6g 6 Qtr[C(R) 2 } 

- 4g 4 C(G)tr[PC(R)} } + g 7 QC(G) [4C(G) - Q] 

= 3r- 1 ^ 3 y ifem y ifen P ri m C( J R)^ + 6r" Vtr [PC(R) 2 } 

+ r-ytr [P 2 C(R)] - 6r -1 Q0 7 tr [C(R> 2 ] - 4r~y C(G)tr [PC(R)] 

+ g 7 QC(G)[AC(G)-Q}. 

(2.25) 
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In the general case, no-one has explicitly computed (3g 3 ^ DRED ; there does, however exist 
a calculation in the special case of a theory without any chiral superfields by Avdeev and 
Tarasov fTTH . For this special case it is easy to see from Eq. ( |2.25| ) that we obtain 

(16tv 2 ) 3 (3^ dred = -21g 7 C(G) s (2.26) 



which precisely agrees with the result of Ref. |T7|] ■ Note that from Eq. ( |2.8|) we have that 



(16n 2 f(3^ NSVZ = -12g 7 C(G) 3 . (2.27) 

This difference was first remarked upon in Ref. Q • The fact that we successfully reproduce 
Eq. (|2.26| ) is an excellent check of both our abelian calculation and our coupling constant 



redefinition. It is intriguing to note that this redefinition, as defined in Eq. fl2.24Q , can be 
written: 

5g = -\pf\ (2.28) 

This of course suggests the possibility of generalising Sg to all orders, and hence deriving 
the all-orders j3^ DRED . We have been unable to provide such a construction, however. 
In the next section we proceed to four loops to test whether the simplicity of our result 
Eq. (|2.28| ) is sustained. 



3. The four-loop calculation 

In the last section we were able to show that the coupling constant redefinition Sg 
relating p^ DRED ^ Q p(3)NSVZ uri jq Ue iy determined (up to an overall constant) without 
any further calculation. Unfortunately this property does not extend to four loops, as is 
easily seen as follows. From Eq. (|2.25|) let us rewrite p^ DRED j n ^ e form 

p{Z)DRED = _ 4Ai + + Aa (31) 



where 



(16tt 2 ) 3 Ai = g 5 C(G) [r" 1 tr[PC(i?)] - g 2 QC(G)} (3.2a) 
(167r 2 ) 3 A 2 = r -1 tr [g 3 S 4 C(R) - 2g 7 QC(R) 2 + 2g°PC(R) 2 ] (3.26) 
(167r 2 ) 3 A 3 = g 3 r-hv[P 2 C(R)] - g 7 Q 2 C(G). (3.2c) 
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The corresponding formula for f3 g ^ NSVZ is 

pf) NSVZ = -4A 1+ 2A 2 . (3.3) 

The purpose of these decompositions is that each of the A, represents a candidate for 
Sg satisfying the requirements we formulated in the previous section; and, indeed, the A* 
are the only such candidates. In particular in the N = 2 case we have Ai = A2 = A3 = 0. 
Therefore we may anticipate that to relate f3 g ^ DRED to (3 g 4 ^ NSVZ we may need to make 
a redefinition of the form Sg = Y2 a i where the on are as yet undetermined coefficients. 
We will also, of course, need to take into account the effect on the four-loop /3-functions of 
the 0(g 5 ) redefinition discussed in the last section. It is immediately clear that we shall 
be unable to determine the coefficient ac\ from an abelian calculation, since Ai vanishes in 
the abelian case. Using Eq. ( |2.18| ), we obtain the leading-order change in f3 g due to each 
of the A^ as follows: 



16tt 2 ) 4 ^(Ai) = 2g 5 C(G)r- 1 jtr [S 4 C(R)] + tr[P 2 C(R)} + 2ghi[PC(R) 2 ] 
+ g 2 Qtv[PC(R)] - 2g 4 Qtv[C(R) 2 ]} - 4g 9 Q 2 C(G) 2 



(3.4) 



(167r 2 ) 4 ,5/3 s (A 2 ) = 2( ? 3 r- 1 tr [(S 7 + 2S 8 + Y*S±Y) C(R) + S 5 P] 
+ 4g 5 r-Hr [{S 4 C(R) + S 5 + P 2 C(R)} C(R)] 
+ 4Qg 7 r-hr [PC(R) 2 - SiC(R)] + 8^ 7 r" 1 tr[PC( J R) ; 
- 8Qg 9 r-hr [C(R) 3 + QC(R) 2 ] 



(3.5) 



(167r 2 ) 4 5/L(A 3 ) = 4g 3 r-Hr \P 3 C(R) + 2g 2 P 2 C{R) 2 - 2g 4 QPC(R) 2 + S 5 P] 

(3.6) 

-4«7 9 Q 3 C(G). 

It is straightforward to verify that each of the d/3 g (Ai) vanishes for N = 2 super- 
symmetry; manifestly they also vanish in a one-loop finite theory. We must be careful 
about the logic of our procedure here, because we know that <y(3)DRED ^ oes no ^ vanish 
in a one-loop finite theory ||18|| , so there is no reason to expect that (3 g 4 ^ DRED will either 
(and, indeed, it does not). There is, however, a redefinition Y — > Y + 5Y corresponding 
to a change to a renormalisation scheme (DRED', say) such that ry(3)DRED ^ oes vanish 

{A)DRED' 



in the one-loop finite case|16[]. In this scheme (5 g will also vanish, relying on the 
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theorem^ that in an n-loop finite theory (3 g n+1 ^ = 0. Suppose that we have found a Sg 
which transforms f3 g ^ DRED to the NSVZ form f3 g ^ NSVZ (i.e. consistent with Eq. (|2.7|) ). 
(We must take into account here the fact that y^NSVZ ^jff ers f rom ^) D R ED j owing 
to the redefinition of g which takes us from DRED to the NSVZ form at the three-loop 
level-see later.) Applying the above SY to j( 3 ) NSVZ and (3^ NSVZ yields a new y( 3 ) NSVZ 
and f3g 4 ^ NSVZ which are also of the NSVZ form, recalling that a redefinition of Y does not 
affect whether or not (3 g satisfies Eq. (|2.7| ). Our redefinition Sg also transforms (3 g ^ DRED 
into (3 g ^ >NSVZ , which both vanish in the one- loop finite case, and hence we conclude that 
Sg must itself vanish in the one-loop finite case. 
From Eq. fl2.7|) , we have that 



(16ir 2 ) 4 P g 4)NSVZ = 8g 9 QC(G) 3 - S/^Gfr^ie^tr [^C(R) 

- Ag b C{G)r- x (l^ 2 ) 2 tT [ 7 (2) C(#)1 (3.7) 

- 2A" 1 (16vr 2 ) 3 tr U^ NSVZ C(R) 

Now in this equation, y^ 1 ' and 7 < - 2 - ) are unambiguous, being defined in Eq. ( |2.3| ) and ( |2.6fc| ) 
respectively, but we must be careful about y( 3 \ 

From Ref. |16| we have the result for ^{3)dred . 

(16vr 2 ) 3 7 ( 3 ) MEB = K {g 6 [12C(R)C(G) 2 - 2C(R) 2 C(G) - 10C(R) 3 - 4C(R)A(R)] 
+ g 4 [4C(R)S l - C(G)S 1 + S 2 - 5S 3 ] + g 2 Y*S x Y + \M 
+ g 2 [C(R)S 4 - 2S 5 - S 6 ] - g 4 [PC(R)C(G) + 5PC(R) 2 ] 
+ 4g 6 QC(G)C(R)} + 2Y* S 4 Y - \S 7 - S s + g 2 [4C(R)S 4 + 4S 5 ] 
+ g 4 [8C(R) 2 P - 2QC(R)P - AQS X - lOr^Tr \PC{R)] C(R)] 
+ g 6 [2Q 2 C(R) - 8C(R) 2 Q + 10QC(R)C(G)] 

(3.8) 

where k = 6£(3). Group theoretic factors are defined in Appendix A. 

Now in the previous section we constructed a 5g that related ^> DRED i Q p^ NSVZ ; 
Eq. (|2.24| ). As we mentioned earlier, this redefinition affects 7 too, so that 

^NSVZ = 7 {3)DRED _ ^SgC^R)^ 2 )- 1 

= yW DRED - 2(16tt 2 )-V [r-Hr [PC(R)} - g 2 QC(G)] C(R). ^ 
This completes the definition of (3 g 4 ' >NSVZ . We anticipate that it will be related to 

^(4) DRED ag follows: 
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(3 (^RED =p (,)NSV Z + J2a i 5f3 g (A i )+n (3.10) 

i 

where O is the change in j3 g ^ due to the redefinition in Eq. ( |2.24j ), 



1 tJY ' 1 t)Y* a c)g\ " " 6g y 

j(5) 



5g-8g.^f). (3.11) 



(There are, of course, terms of 0((<5g) 2 ), but these are 0(g 11 ) and hence affect (5 g 
Substituting for f3y and j3 g 2 ^ we obtain: 

(16tt 2 ) 4 = g 3 r- 1 {2Qghv[PC(R) 2 + 2g 2 C(R) 3 ] - 2g 2 tv[P 2 C(R) 2 + 2g 2 PC(R) 3 ] 

+ 2# 4 Qtr [SiCiR)] - tr[Y*S A YC(R)} 

- tr [2g 2 (S 5 + S A C(R)) C(R) + S 5 P] }. 

(3.12) 

Given an explicit calculation of (3^ DRED ', we would be able to test the validity of the 
construction Eq. fl3.10|) ; and if it proved to be valid, determine the o^. This calculation is 
beyond our strength; we have, however, performed a partial calculation consisting of the 
contributions to f3g ED in the abelian theory of 0(g 3 Y 6 ). From Eqs. ( |3.4|) - (|3.6| ) it is 
clear that this will enable us both to test whether our construction works with regard to 
such terms and also to determine a 2 an d «3, thus fixing (3 g ^ DRED apart from a single 
unknown parameter (a±). The calculations were performed using the methods explained 
earlier in the three-loop case. The relevant diagrams are shown in Fig. 3, and the results 
are given in Table 2. As before, the momentum integrals can be expressed in terms of 
a convenient basis, depicted in Fig. 4. The simple poles for each of these (subtracted) 
diagrams are given by 

-^simple = 777-2^47 ( o ~ ^( 3 ) ) ' ^simple = 2 777^2717 ( _1 + ^( 3 ) 



(167r 2 ) 4 e V2 M V p (167r 2 ) 4 e 

1/5 \ 

-^simple 777 n\4 ( 7 C(3) J j -^simple 



(167r 2 ) 4 e V 6 ^ 7' bimpie 2(167r 2 ) 4 e' 

T - 2 1 K - 5 1 

Simple 3 (167r2) 4 e > g (16^2)4, 

I 1 _ 11 1 

^simple- 2(167r 2)4 e ' ^simple - g^g^,' 

II 11 

> simple — 7 7T7 TTT" ' -'simple 



6(167r 2 ) 4 e' " mpie 6(167r 2 ) 4 e 

^simple 3(167r 2 ) 4 e > ^simple- 2(^2)4^ 



(3.13) 
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It is interesting to note that £(3) dependence only appears in those diagrams which con- 
tain a "figure-of-eight" sub-diagram. This could perhaps be explained by invoking recent 
proposals which relate the appearance of transcendental numbers in Feynman integrals to 
knot theory [|l~9f . On adding the results for all the diagrams, we find 

{\Qn 2 fd^ DRED = y 3 r-Hr [(2P 3 - 2S 8 - 19Y*S 4 Y - S 7 ) C(R)] 
- \Kg 3 r-Hi[MC(R)] - y 3 r-Hi[S 5 P] 

+ ■ ■ ■ (terms involving C(G), and of order g 5 Y , g Y 2 and g 9 .) 

(3.14) 

The coefficient of the tr[MC(R)] term was derived not by direct calculation, but by ex- 
ploiting the fact that, as we observed earlier, the coupling constant redefinition which 
makes 7^ 3 - ) vanish in the one-loop finite case should also make (3g^ vanish. It was shown 
in Ref. [jTIJ that the redefinition 

(167r 2 ) 2 SY ljk = lKY Um Y jpq Y krs Y lpr Y mqs + 0(g 2 ) (3.15) 

makes 7 < - 3 - ) vanish in the one-loop finite case. Since a change 5Y induces a leading-order 
change in f} g of the form 



6/3 B 



d d 
SY.— + SY*.- 



dY dY 
(16ir 2 )~ A r~ 1 g 3 C(R) i j(Y jkl 6Yiki + SY^ kl Y tkl 



8^ 

9 (3.16) 



we can deduce the coefficient of tr[MC(R)] to be as given. This indirect deduction was 
necessitated because the direct calculation of this coefficient proved very involved. One 
might feel slightly uneasy because it is not clear that the theorem that flg n+1 ^ vanishes in an 
n-loop finite theory should be valid in an arbitrary renormalisation scheme. To allay these 
doubts, we have explicitly computed the coefficients of tr[C (R) 2 A(R)] and tr[Y* SiYC(R)} 
m $ )DRED in the one-loop finite case, and checked that they agree with those obtained 
by the same indirect argument. 

Now comparing the result Eq. ( |3.14| ) with Eqs. (|3.4|) - ([3.7|) and ( |3.12| ), we see that we 
require 

« 2 = -§, « 3 = \. (3.17) 

We note that it is non-trivial that a solution exists at all for «2 and 0:3, since we need 
to reproduce six terms in Eq. (|3.14| ) and we have only two free parameters «2 and as to 



adjust. Note, however, that the coefficient of tr[MC(R)] automatically satisfies Eq. (|3.7|), 
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irrespective of the values of a 2 and a^; this is essentially guaranteed, as we know that we 
can redefine Y so that the M terms vanish in both 7^ and f3^ (hence trivially satisfying 
Eq. (|3.7|) as far as these terms are concerned), and we also know, as mentioned earlier, 
that redefinitions of Y have no effect on whether Eq. ( |3.7f ) is satisfied. 

As we have indicated, several miracles were required to facilitate our construction. It 
is perhaps disappointing, however, that it remains unclear how the redefinition we have 
found generalises to higher orders. From Eq. ( |2.28|) , it was tempting to conjecture that 
the transformation Sg = c^-i^i might have been proportional to p^ DRED or p^ NSVZ 



but it is easy to see from Eqs. ( p7T|) and (|3.3|) that this doesn't work. 
The final complete result is 

(16tt 2 ) 4 /3^ dred = g 3 r-Hv [{-\kM - fY*S 4 Y + §P 3 - ±S 7 - §S 8 } C(R) - |5 5 P] 
+ 5 r _1 tr ^-2kY*S 1 Y - (2k + f)C(R)S 4 + (4k - f)S 5 + 2kS 6 
- ^P 2 C(R) + (2c*i + A)C(G)S 4 + 2a 1 C(G)P 2 }c(R) 

+ g 7 r-hr ^-8kC(R)S 1 + 2kC(G)S 1 - 2kS 2 + WkS 3 
+ 2Ar-hi[PC(R)}C(R) + (10k - f)PC(R) 2 + 2QPC(R) + fQS 1 
+ (2k + 4ai + 8)C(G)PC(R) + 2 ai QC(G)P - 8C(G) 2 p}c(R) 



+ g 9 r 1 tr 



f QC(Rf + ±Q 2 C(R) 2 - (8k + 4ai + 32)QC(G)C(R)' 



- 2AkC(GYC(RY + 4kC(G)C(RY + 20kC(RY + 8kA(R)C(R) 
+ g 9 [8QC(Gf -Aa 1 Q 2 C(G) 2 - |Q 3 C(G)] 

(3.18) 

We see that in a one-loop finite theory (P = Q = 0), we have j3^ DRED 7^ 0; but, 
as noted earlier, it can be transformed to zero by means of the redefinition Y — > Y + 8Y, 
where SY is the transformation that makes ^(3)dred y^jg^ j n guch a theory. We gave 
the 0(Y 5 ) term in this 5Y in Eq. ( |3.15| ); the full expression is[TE 



(16iv 2 ) 3 dY ljk = K{\Y ilm Y m Y krs Y lpr Y mqs + \g 2 m Y jk)m 

+ g 4 [C(R) (l rn C(Ry n Y^ mn - lY n ^ k C(R)^ m C(R) m n - A(R)Y ljk 
- \C(G)C(R) {l rn Y jk)m +3C(G) 2 Y ljk }}. 

(3.19) 

This redefinition defines the DRED' scheme which we introduced previously; correspond- 
ingly, the same redefinition applied to (3g 4 ^ NSVZ ', as given in Eq. (|3.7|) , will define the 



NSVZ' scheme for which (given P = Q = 0) f3 { g 4)NSVZ ' = likewise (as can readily be 
checked using Eq. ( |3.16| )). 
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4. Conclusions 

We have explicitly constructed the coupling constant redefinitions that relate j3^ svz 
to f3g ED up to and including four loops, except for one undetermined parameter. The fact 
that this construction was possible demonstrates that the renormalisation scheme in which 
the NSVZ form is valid is perturbatively related to the conventional DRED scheme. As a 
by-product of our investigations we have obtained p^ DRED anc [ f3^ DRED for a general 
non-abelian theory, except for the dependence of f3^ DRED on the same undetermined 
parameter. 

The fixed points of (3^ svz , which satisfy the equation 

2r _1 tr [-yC(R)] = Q (4.1) 

are significant in dual gauge theories. Now the existence of a fixed point is preserved under 
a coupling constant redefinition g — > g'(g), as long as the function g'(g) is differentiable 
at the fixed point. Our demonstration that the NSVZ and DRED schemes are perturba- 
tively related therefore suggests that there is a fixed point of [3 ERED corresponding to any 
fixed point of (3^ svz . It would have been interesting had we been able to construct the 
DRED «-> NSVZ redefinition to all orders; unfortunately, however, to the extent that 
we have pursued the perturbative form, there is no indication of an all-orders result. It is 
tantalising, in this regard, that in the development of background-field covariant superfield 



Feynman rules |2(J it appears that beyond one loop, loops of e-scalars play a crucial role 
in the relevant counter-terms. Such a loop provides a factor of e, so that the simple pole in 
e thus depends on what would have been the double pole had the e-scalars been ordinary. 
Given that all higher order poles in e are determined in terms of the simple pole by the 
't Hooft consistency conditions [^T|, one might have hoped to proceed to an all-orders 



construction of f3 ERED . In this endeavour we have not been successful. 

Nevertheless, we feel that the exercise has been worthwhile. We have demonstrated be- 
yond all reasonable doubt that there does exist a scheme in which the NSVZ /3-function is 
valid. Our result for (3^ DRED , in conjunction with the result for ^(^dred f rom j^ e f j 



is in any event of interest phenomenologically ||22| |, especially, perhaps, in post-post-modern 



theories with additional matter content (see for example |[23| ). 
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Appendix A. Notation and Conventions 

Here we give some details concerning our notation and various useful formulae. With 
our conventions, the -D-algebra in Minkowski space is 

{A*, A*} = faaa 9 ^ 

_ _ (A.l) 
{D a ,Dp} = {D a ,D } = O, 

where a M = (J, a 1 ), a 1 being the Pauli matrices. It is then easy to verify that the projection 
operators 

n = D 2 D 2 + D 2 D 2 , 

(A.2) 

III = -2D a D 2 D a 

2 

satisfy 

ni+n = -<9 2 (A.3) 

as stated in Section 2. We use standard superfield Feynman rules, as described, for ex- 
ample, in Ref. J24| , with minor changes of normalisation consequent upon our conventions 
in Eq. (|A.l] ); the only unusual feature being the way we use transversality of the vec- 
tor self-energy to simplify the D-algebra, as described in Section 2. A useful identity in 
manipulation of supergraphs (which follows immediately from Eqs. (|A.2| ), ( |A.3| )) is 

D 2 D 2 D 2 = -d 2 D 2 , D 2 D 2 D 2 = ~d 2 D 2 . (A.4) 

Gauge invariance of the superpotential means that the Yukawa couplings Y l ^ k satisfy the 
equation 

Y m ^R k 2m = (A.5) 
From this identity it is straightforward to prove that: 

Y imn B? Am R k An = \ [Y^ k C(RY m - Y^ k C(Ry m - Y^ m C(R) k m ] . (A.6) 

Eq. ( |A.5| ) is rather similar to momentum conservation at a three-point vertex; correspond- 
ingly, Eq. (|A.6|) is analogous to the simple identity 

P-q=l[(p + q) 2 -p 2 -q 2 } (A.7) 
17 



Eq. ( |A.6|) is very useful in dealing with the group theoretic factors. It is always possible, in 
the calculations we have performed, to push Ra's around so as to produce the quadratic 
Casimir C(R) (though there is no reason to expect this property to persist to arbitrary 
orders). In this process, much labour involving dummy indices is avoided by the adoption 
of a diagrammatic notation where Y or Y* are represented by vertices, index contractions 
5 % j by propagators, and Ra by "mass insertions" . 
We use the following definitions: 

= Y imn C{RY m Y 3vn (A.80) 

(Y*S 1 YY J = Y^S^mYjpn (A.86) 

S^ = Y irnn C{Ry m C{R) q n Y 3pq (A.8c) 

S^ = Y* mn (C(R) 2 ) p m Y Jpn (A.8d) 

= Y P p m Yjp n (A.8e) 

(Y*S A Yf 3 = Y^ n S 4 p m Y Jpn . (A.8/) 

Sis = Y imn C{RY m P\Y 3nq (A.8g) 

Sis = Y imn C{RY m P\Ys Pq {k.d,h) 

qi . _ Y imn PP pi V. ( A 

Djj — I r m r n I jpq \rt..Ol) 

Sis = Y* mn (P 2 y m Y Jpn (A.8j) 
A(R) = Y,C(R a )T(R ct ) (AM) 

a 

= Y M Y kmn Y lrs Y pmr Y qns Y Jpq (AM) 
In Eq. (|A.8Aj ) the sum over a is a sum over irreducible representations. 
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Diagram 


Integrals 


Group factor 


a 







b 


C 


2X 4 + 4X 3 - 2Xi 


c 


A 


4Xi 


d 


B 


2X 4 


e 


B 


4X 4 


f 


2A + 2B 


-I6X3 


g 


A 


-8X3 


h 


A 


I6X3 


i 


B 


I6X3 


j 


B 


I6X3 


k 


-2A+\C 


8g 6 Qtr[C(R) 2 } 


1 


-\A+\B 


16g 6 Qtr[C(R) 2 ] 


m 


-\A+\B 


-32g 6 Qtr[C(R) 2 } 


n 


B 


8g 6 Qtr[C(R) 2 ] 


o 


-B-2D 


8X 3 + 2X 4 + 8gHr[C(R) s ] 


P 







q 


A 


-8X 2 + 24X 3 + 48# 6 tr[C( J R) 3 ] 


r 


2A + 2D 


4X 2 - I2X3 - 24# 6 tr[C( J R) 3 ] 


s 


2A-2D-E 


4X 2 - I2X3 - 24# 6 tr[C(i?) 3 ] 


t 


A 


8X3 - 8X 2 + 16p 6 tr[C(i2) 3 ] 


u 







V 


E 


4X 2 - 4X 3 - 8gHr[C(R) 3 ] 


w 


2A + B + 2D 


8X3 + 16^ 6 tr[C(i?) 3 ] 


X 







y 


AA-C -E 


8X3 + 16/tr[C(i?) 3 ] 


z 


A 


-I6X3 - 32g 6 tr[C(R) 3 ] 



Table 1: Three-loop contributions to < VV > po i e for the diagrams in Figure 1. Each 
contribution is obtained by multiplying the simple pole from the momentum integral in the 
first column by the group theory factor in the second column, and by r~ l . 



Diagram 


Integrals 


Group factor 


aa 


A 


16X 3 + 32gHr[C(R) 3 ] 


bb 


2A + 2B 


-32qHr\C(R) 3 ] 


cc 


B 


16qHr\C(R) s ] 


dd 


B 


16g 6 trfC(i?) 3 l 

,) L \ J J 


ee 


A 


32 9 6 tr[C( J R) 3 l 


ff 


B + E 


-32qHr\C(R) 3 } 


Kg 


E 


16g 6 trfC(i?) 3 l 


hh 


4A + 3B - 2D 


8g 6 tr[C(i2) 3 l 


ii 


2B - 8D + 2E 


8q 6 tr\C(R) 3 } 


ii 

jj 


AB + C + 2D + E 


16g 6 trfC(i2) 3 l 


kk 


B + E 


16g 6 trfC(i?) 3 l 

l7 L V / J 


11 


B + E 


-32g 6 tr[C(i2) 3 l 


mm 


B 


32o 6 tr[C( J R) 3 l 


nn 


A 


32g 6 tr[C( J R) 3 l 


ii 


2A + 2B 


-32g 6 tr[C(i2) 3 l 


QQ 
mm 


B 


-16g 6 tr[C(.R) 3 l 


rr 


2A - 2D 


-16g 6 tr[C(.R) 3 l 


ss 


A 


-32g 6 tr[CCR) 3 l 


tt 


A 


32g 6 trfC( J R) 3 l 

L \ / J 


uu 


B 


16g 6 tr[C(i?) 3 ] 


vv 







WW 







XX 








Ta&te i (continued) 



Diagram 


Integrals 


Group Factor 


a 


4P-2R + 2J 


tr[S 7 C(R) - 2S 5 P] 


b 


F + 2N 


4tr[P 3 C{R)} 


c 


J + 2Q - L 


2tr[S 7 C(R)] 


d 


J + 2P 


8tr[S 5 P] - Atr[S 7 C(R)] 


e 


-I 


Atr[Y*S 4 YC(R)] 


f 


-M 


2tr [-Y* S 4 YC(R) + S 5 P] 


o 







h 







i 







i 

j 







k 


—L 


-4tr\S 7 C(R)} 

L / \ /J 


1 


-G 


2trf(P 3 - S 8 )C(R)} 


m 


-K 


2tr\S 7 C(R)} 

L * V / J 


n 


-H 


4tr[S 8 C(i2)l 

L o V /J 


P 


-J 


8tr[5 5j P] 


q 


-J 


4tr[5 5 P] 


r 


-F 


4tr[P 3 C(i?)] 


s 


-F 


4tr[P 3 C(i?)] 



Table 2: Four-loop contributions to < VV > po i e for the diagrams in Figure 2. Each 
contribution is obtained by multiplying the simple pole from the momentum integral in the 
first column by the group theory factor in the second column, and by r~ 1 g 2 . Note that only 
contributions including terms of 0(Y 6 ) have been retained in the group theory factors. 




(XX) 

Fig. 1. (Continued) 




Fig. 2. Three-loop momentum integrals 




Fig. 4. Four-loop momentum integrals 



